Sequential Quantum Cloning 
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Not all unitary operations upon a set of qubits can be implemented by sequential interactions 
between each qubit and an ancillary system. We analyze the specific case of sequential quantum 
cloning, 1 — * M, and prove that the minimal dimension D of the ancilla grows linearly with the 
number of clones M. In particular, we obtain D = 2M for symmetric universal quantum cloning 
and D = M + 1 for symmetric phase-covariant cloning. Furthermore, we provide a recipe for the 
required ancilla-qubit interactions in each step of the sequential procedure for both cases. 

PACS numbers: 03.67.-a, 03.67.Mn, 03.67.Dd 



Multipartite entangled states stand up as the most 
versatile and powerful tool to perform information- 
processing protocols in Quantum Information Science [l[ . 
They arise as an invaluable resource in tasks such as 
quantum computation 0, 01, quantum state teleporta- 
tion 01, quantum communication 0| and dense coding 
As a result, the controllable generation of these 
states becomes a crucial issue in the quest for quantum- 
informational proposals. However, the generation of mul- 
tipartite entangled states through single global unitary 
operations is, in general, an extremely difficult exper- 
imental task. In this sense, the sequential generation 
studied by Schon et al. Q , where at each step one qubit 
is allowed to interact with an ancilla, appears as the most 
promising avenue. The essence of this sequential scheme 
is the successive interaction of each qubit initialized in 
the standard state |0) with an ancilla of a suitable di- 
mension D to generate the desired multiqubit state. In 
the last step, the qubit-ancilla interaction is chosen so as 
to decouple the final multiqubit entangled state from the 
auxiliary D-dimensional system, yielding Q 
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Here, the VjjjJ are D— dimensional matrices arising from 
the isometries (unitaries) Vjfc] : \)a ® (|0)) — > \)a ® §B k , 
with f)yi = C D and f)B k = C 2 being the Hubert spaces 
for the ancilla and the fcth qubit, respectively, and where 
\tpi) and \(Pf) denote the initial and final states of the 
ancilla, respectively. The state ([1]) is, indeed, a Matrix 
Product State (MPS) (cf. e.g. [8|and references therein), 
already present in spin chains [ 9|, classical simulations 
of quantum entangled systems |10fl a nd density-matrix 
renormalization group techniques [111 ]. Moreover, it was 
proven that any multiqubit MPS can be sequentially gen- 
erated using the recipe of Ref. Notice that in this 
formalism, the mutual qubit-ancilla interaction in each 
step k completely determines the matrices V?£ , ik = 0, 1, 



whereas we enjoy some freedom to build such an inter- 
action from a known Vro . This freedom stems from the 
fact that in the proposed scheme only the initial state |0) 
for each qubit is relevant. 

In this letter, we consider the possibility of implement- 
ing quantum cloning based on a sequential protocol with 
the help of an ancillary system. This problem is cer- 
tainly far from being an application of Ref. 0] , given that 
the initial and final states are unknown. In this sense, 
any proposed strategy will be closer to the open problem 
of which global unitary operations (certainly not all of 
them) can be implemented through a sequential proce- 
dure. Despite the fundamental no-cloning theorem 
stating the impossibility to exactly clone an unknown 
quantum state, there exists several cloning techniques 
with a given optimal fidelity [l3| . These procedures dif- 
fer either from the initial set of states to be cloned or from 
symmetry considerations. In general, an optimality con- 
dition of the cloning procedure is obtained via the max- 
imization of the fidelity between the original qubit and 
each final clone state. We will show how to perform se- 
quentially both the universal symmetric [3, Il5| and the 
economical phase-covariant symmetric quantum cloning 
[lil , [ij} from one qubit to M clones. In the first case, a 
global unitary evolution transforms any input state 
in a set of M clones whose individual reduced states p out 
carry maximal fidelity with respect to \tp): fi m — 2A 3 I A ^ 1 ■ 
This cloning procedure is fully described by the evolution 

|V)®|-B) \GM M (ip)) = 

M-l 

= X>J'K M - 3)1>,M> X )s®\{M- j - l)V'*,^* X )s, 

3=0 

(2) 

where |GMm(V0) stands for the state produced by the 
Gisin-Massar cloning procedure [H|], that results in M 
optimal clones of from the initial blank state \B), 



2 



o-j = \J m(m+{) , and \(M - fiipjip^s denotes the nor- 
malized completely symmetric state with (M — j) qubits 
in state ip and j qubits in state i/j . Notice the presence 
of M — 1 additional so-called anticlones. They are nec- 
essary in order to perform this cloning procedure with 
the optimal fidelity. The anticlone state tp* refers to 
the fact that they transform under rotations as the com- 
plex conjugate representation. For concreteness sake we 
have chosen = cos#/2|l) + e~^sin#/2|0) in coinci- 
dence with the seminal paper by Buzek and Hillery [Tij ]. 
whereas |V>) = cos6»/2|0) + sin0/2|l). In the second 
case, motivated by quantum cryptoanalysis, the goal is to 
clone only those states belonging to the equatorial plane 
of the Bloch sphere, i.e. those such that 6 = ir/2. Fur- 
thermore, we have only focused upon the cases where 
no anticlones are needed (hence the term economical). 
Under this assumption, imposing the purity of the joint 
state, the number of clones M must be odd [l(|. The 
cloning evolution is now given by 

\tP) ® \B) -h. -L [\(k + 1)0, kl) s + e l0 |fcO, (k + l)l) s ] , 

(3) 

where k = (M — l)/2 and where we have followed the 
same convention as above. 

The basic idea is to express the final states ((2|) and 
([3|) in its MPS form, as given in Ref. [h}, by performing 
n — 1 sequential Schmidt decompositions 

l*>= E I^L 1 l>A[l] ai lA 2 >---A[n-l] Q „_ 1 |^Ll I )> 

ai ...a n — ± 

and then writing the unnormalized Schmidt states 
in the computational basis for the corresponding 
qubit \^_ iai ) = Ein}%. iai \il). Then, |$) = 

Ei 1 ...i jv c u-wh'i---*jv) ! with 

c il ...i„= E r[i]^ i A[i] Q1 r[2]^ iQ2 A[2] Q2 ...r[n]j rii _ i . 

ai ...a n — i 

(4) 

We identify the matrices Vj^ by matching indices in ex- 
pressions (HJ) and ([4]). The indices aij run from 1 to x, 
where x = m ax-p{x-p}, Xv denoting the rank of the re- 
duced density matrix p-p for the bipartite partition V of 
the composite system [lOj. 

In order to employ the sequential ancilla-qubit device 
as a quantum cloning machine we will firstly elucidate 
the minimal dimension required for the ancilla. To clone 
an arbitrary input qubit state = a|0) + /3|1), we ex- 
ploit linearity and determine the minimal dimension Do,i 
of the ancillas to perform the cloning for the states |0) 
and |1) and then combine both results in a single ancilla 
of minimal dimension D to be determined. Let us focus 
upon the symmetric universal cloning of |0). To deter- 
mine the minimal dimension Dq of the ancilla we need to 
compute x-i which can be done without the exact MPS 
expression for the state. 



Let us denote by V = A\B the partition into two 
subsystems, one with the first A qubits, the other 
with the following B qubits, and Ca\b the correspond- 
ing coefficient matrix. For definiteness, CxibCVO = 
[ci 1 ...i A ,i A+1 ...i A+B \i where i\...iA is treated as the row 
index, whereas ia+i ■ ■ ■ iA+B is treated as the column in- 
dex, and Ci 1 ...i A ,i A+1 ...i A+B denote the coefficients of state 
| t/j). Now, the Gisin-Massar state cloned from |0) can be 
written as 

|GM M (0)) = 

Af-l 

S M ® 5m-iE a i\( M - ^J 1 ) ®\{M-j- 1)1, jO), 

(5) 

where Sa is the normalized symmetrizing operator for A 
qubits, so that Sm ® Sm-i is an invertible local oper- 
ator for the partition M\M — 1. Due to the orthonor- 
malities among the states on the rhs, their Cm\m-i can 
only have M different rows whereas the rest are all null, 
hence r(C M | M _ 1 ) = M. As Sm g) 5a/_i amounts to 
local changes of basis within both partitions only, they 
cannot change the rank of the density matrix Pm\m-i, 
so that the rank of the coefficient matrix of @ is also 
M. We now consider the partition k\2M — k — 1, where 
k = 1, . . . M — 2. The matrices C^M-k-i ar e obtained 
from the Cm\m-i by adjoining rows and columns to make 
them longer, but - as that there are only M different 
rows in Cm\m~i, the rest being all null - this reordering 
procedure cannot increase the former rank. Finally, 

r(C kl2M -i- k ) < r(<?M|M-i) = M. (6) 

From the results above, it follows that x = M, i.e. 
that the minimal dimension Dq to clone the |0) state 
is Dq = M, namely the number of clones to produce. 
Repeating the same argument for the initial state |1) we 
also conclude that the minimal dimension of the ancilla to 
clone the |1) state is D\ = M, as expected. Now we must 
combine both results to find D for an arbitrary unknown 
state = a\0) + f3\l). It is a wrong guessing to think 
that it should also be D = M and, consequently, a differ- 
ent scheme must be given. The MPS expression of @ for 
the original state |0) determines the Z)-dimcnsional ma- 
trices Vqh.] i whereas the corresponding MPS expression 
for the original state |1) determines V\uai 

\gm m (o)) = ■ ■ ■ Kk ■ ■ ■ 

ii..4„=0,l 

\GM M (1)) = £ \V^ n] . . . V*fa \0} d \h . . . i n ). 

ii...i n =0,l 

(7) 

Here, |<^') and \<fp ) can be calculated explicitly and 
will play an important role below. 
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We propose now to double the dimension of the ancilla, 
C D — > C 2 ® C"°, in order to implement a deterministic 
protocol of sequential quantum cloning. 



Protocol 1. i. Encode the unknown state 
initial ancilla state \tpi) = \ip) ® \Q)d- 



in the 



ii. Allow each qubit k to interact with the ancilla ac- 
cording to the 2D -dimensional isometries Vj]S — 

|0><0|®^ ]+ |1><1|®^ 



rtk 
L[fe]- 



Perform a generalized Hadamard transformation 
upon the ancilla 



i (o)\ 1 



|o> ® |^ 0) ) 



1^ 



(8) 



Note that the choice C D — > C 2 ® C D (based on ped- 
agogical reasons) could be changed, equivalently, to 
C D — > C 2D . In this way, Eq. {5p would not display 
entangled states but simple linear superpositions. 

iv. Perform a measurement upon the ancilla in the lo- 
cal basis {|0) ® |1) ® l^ 1 ')}. 



v. If the result is \Q) ®\(fp) (which happens with prob- 
ability 1/2), the qubits are already in the desired 
state; if the result is |1) ® \cpp ) (probability 1/2), 
perform a local ix-phase gate upon each qubit, then 
they will end up in the desired state. 

Proof. After the first two steps, the joint state of the 
ancilla and the qubits is a (jO) ® \tp^)^ ® \GM M (0)) + 

(3 (Jl) ® 1^)) |GM M (1)), where originally \ip) = a\0) + 
/3\1). After the Hadamard gate in (iii), this state becomes 

-J= (|0) ® |^ 0) )) ® [a|GAfjif(0)) +/3|GM M (1)>] + 
+ i= (|1) ® l^)) ® [«|GM M (0)) -/3|GM m (1)>] . 

The remaining steps follow immediately from this expre- 
sion and from linearity [la ]. □ 



Notice that despite the measurement process in step 
(iv), the desired state is obtained with probability 1, 
while the fidelity of each clone is optimal, F^m = 2 ^j" 1 , 
as in Ref. flEj . In summary, the minimal dimension D 
of the ancilla for cloning M qubits is D = 2 x M, i.e., 
it grows linearly with the number of clones even if their 
Hilbcrt space grows exponentially (2 M ). 



k = 1 



M!i]]<,= 



5ijC(2 



1) 1 < t, j < 2 
otherwise 



8 l , 3 _ J C(2-i,i~-l) l<i,j<2 
-j^Sij otherwise 



o[»]J 



C(n+l-i,i-l) 



1 < i,j < n 



J ^ C(n-i,i-l) 

-j^Sij otherwise 



i 

C(n-j,j) 



i = 1 ; j = n + 1 
7+1 ^"-'^ 2<i<n + l;l<7<n 

'J + 1 C(n-3j-l) — — 

ot her wise 



1 r 



C(Af-i,i-l 



1 < i, j < M 



Oi,j+l 



[A/ + n]J 



o 

73 



1 < i < M - n 
2<j<M-n + l 
i = M - n+ 1; 1 < j < AI 
otherwise 



r. . / A/-n+l-i 
M-n 



1 r 
V2 ° i] 



1 < i,j < M - n 
i = M - n+ 1; 1 < j < AI 
otherwise 



TABLE I: Matrices for the universal symmetric cloning protocol. 



It can be checked straightforwardly that if one had to clone a d-dimensional system, the minimal dimension for 
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the ancilla would be D = d x M, an obvious generaliza- 
tion of the preceding argument. 

For the symmetric phase-covariant cloning, the same 
arguments can be reproduced. For example, the first 
term on the r.h.s. of Eq. ^ can be cast in the form of 
the state in Eq. @ 

|(fc + l)0,M) s = 

k 

= 1 - 3)0, n)s ® \(k-j)l,jO) s , (9) 

where jj ^ for all j, and similarly for the second term. 
Thus for symmetric phase-covariant cloning the minimal 
dimension for the ancilla is D = 2 x (fc + 1) = 2 x = 
M + 1. We see that the dimension of the ancilla D also 
grows linearly with the number of clones, although it is 
now lesser than above. This is a direct consequence of the 
reduction in the set of possible original states to clone. 
For the symmetric universal cloning we give in detail 
in Table |l] the 2D— dimensional matrices Vr^ driving us 
to a concrete sequential scheme, and where C(i,j) = 




— if q > p and 



1 < n < M — 1. Furthermore, we also have = 
where by i we indicate i = £®1 (mod 2). They coin- 
cide also with the ones for the symmetric phase-covariant 
cloning just by doing the substitutions M — > M 2 I ~ 1 and 

It can be readily verified that the minimal dimension for 
the ancilla is 2 x M. When sequentially applying these 
matrices to the initial state \(pj) of the ancilla, one can 
check, as expected, that if we were to stop at the Mth 
step, the M clones would have already been produced 
with the desired properties, although in a highly entan- 
gled state with the ancilla. To arrive at a final uncoupled 
state, the remaining M — 1 anticlones must be operated 
upon by the ancilla. Note the exponential gain achieved 
with this protocol; despite the 2 M -dimcnsional Hilbcrt 
space of the M clones, we just need a 2Af-dimcnsional 
ancilla. This is a consequence of the Matrix-Product de- 
composition of the Gisin-Massar universal cloning state. 
The proposed schemes can be implemented in a variety 
of physical setups: microwave and optical cavity QED, 
circuit QED, trapped ions, and quantum dots, among 
others. As a paradigmatic example, the clone could be 
codified in a photonic state and the ancilla in a D-level 
atom 0] , and the sequential operations carried out by Ra- 
man lasers would produce unitaries associated with the 
isometries appearing in Table I. These and other re- 
quired unitary operations, as local Hadamard gates, are 
standard in most of the above mentioned physical setups, 
making our proposal suitable for future implementation. 

In conclusion, we have shown how to reproduce sequen- 
tially both the symmetric universal and symmetric phase- 



covariant cloning operations. For the universal cloning 
we have proved that the minimal dimension for the an- 
cilla should be D = 2M, where M denotes the number of 
clones, thus showing a linear dependence. The original 
state must be encoded in a 2 D— dimensional state. For 
the phase-covariant case, the required dimension D of 
the ancilla can be reduced to D = M + 1. In both cases, 
the ancilla ends up uncoupled to the qubits. Along simi- 
lar lines, this sequential cloning protocol can be adapted 
to other proposals, such as asymmetric universal quan- 
tum cloning machines or other state-dependent proto- 
cols. This procedure can have notable experimental in- 
terest, since it provides a systematic method to furnish 
any multiqubit state using only sequential two-system 
(qubit-ancilla) operations. 
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